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CHAPTER  1‘ 

INTRODUCTION 

In  this  study  we  will  be  concerned  with  the  development  of  a 
hyperstable  model  reference  adaptive  control  for  a class  of  nonlinear  plants. 
These  results  represent  an  extension  of  the  Landau  [19]  method  for  the  linear 
time  Invariant  plant  to  the  nonlinear  case. 

In  recent  years  there  has  been  much  engineering  Interest  In  adaptive 
control  systems.  One  of  the  reasons  for  this  Interest  Is  that  modeling  of 
systems  often  Introduces  parameter  value  and  structure  errors.  When  these 
errors  are  carried  through  to  the  control  design  stage,  they  often  produce 
unacceptable  controls  for  the  real  system.  Hence,  a controller  which  "adapts" 
to  these  errors  so  as  to  stabilize  the  dynamic  characteristics  of  the  system 
about  the  specified  response  would  be  highly  advantageous.  This  type  of 
control  was  believed  to  have  been  first  developed  by  Whitaker,  Yamron,  and 
Kezer  [1]  In  1958.  Since  then,  their  Initial  work  has  motivated  an  adaptive 
control  theory  encompassing  many  different  adaptation  approaches.  It  has  been 
the  success  of  these  techniques  to  problems  ranging  from  missile  flight 
control  to  biological  systems  to  economics  and  source  management  [2],  that 
motivate  this  study. 

As  Is  suggested  above,  the  adaptive  control  problem  often  arises 
when  exact  modeling  Is  Infeasible.  However,  If  the  main  source  of  error  Is 
of  the  structure  type,  little  can  presently  be  done  except  for  some  special 
cases.  Hence,  the  discussion  Is  limited  to  parameter  adaptation.  The 
following  background  information  is  presented  to  expose  the  underlying 
problem  and  mathematical  treatment  Involved. 
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The  parameter  adaptation  problem  has  many  sources.  As  an  example, 
consider  the  following  typical  control  problem:  Given  a plant  with  some 

inputs  and  outputs,  design  a controller  such  that  the  plant  responds  in  a pre- 
specified manner.  To  solve  this  problem,  one  would  probably  first  try  modeling 
the  plant  process  with  some  mathematical  structure,  such  as  a set  of  differential 
equations,  which  together  with  some  uniquely  defining  parameters,  would  reason- 
ably describe  the  physical  behavior  of  the  plant.  One  would  then  search  for 
a controller  for  this  model  and  then  hope  that  the  same  controller  would  work 
in  the  plant.  However,  this  is  often  an  infeasible  solution  for  the  following 
reason:  Parameter  errors  will  exist  whenever  the  modeling  fails  to  take  large 

parameter  variations  into  account.  These  parameter  variations  often  take  the 
physical  form  of  parts  wearing  out,  manufacturing  tolerances  in  the  machine 
parts,  or  as  physical  variations  due  to  the  different  operating  conditions  of 
the  equipment  (as  is  the  case  of  the  space  shuttle).  The  net  effect  in  some 
of  these  cases  may  be  unacceptable  performance  or  instability  in  the  plant. 

In  fact,  under  some  circumstances  the  initial  parameter  values  may  be 
indeterminable  and  time  varying  so  the  conventional  design  may  be  infeasible. 

Thus  parameter  errors  are  inherent  in  some  control  applications.  The  question 
is  whether  they  can  be  neglected  in  a particular  case,  and  if  not,  what  can 
be  done  about  it . 

In  some  cases,  the  solution  lies  in  parameter  adaptive  control,  which 
for  this  study,  will  be  defined  as  in  Landau  [3] : 

"An  adaptive  system  measures  a certain  index  of  performance  using 
the  inputs,  the  states,  and  the  outputs  of  the  adjustable  system.  From  the 
comparison  of  the  measured  index  of  performance  values  and  a set  of  given  ones. 
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the  adaptation  mechanism  modifies  the  parameters  of  the  adjustable  system  or 
generates  an  auxiliary  input  in  order  to  maintain  the  index  of  performance 
values  close  to  the  set  of  given  ones." 

Thus,  if  such  an  adaptive  controller  could  be  realized,  one  would 
be  able  to  solve  the  above  control  problem  despite  the  existence  of  parameter 
value  uncertainty.  The  case  with  unknown  parameters,  though  actually  more 
difficult  to  treat  mathematically,  would  similarly  be  solved. 

The  methodologies  required  to  analyze  and  synthesize  the  adaptive 
controller  falls  into  one  of  two  categories:  the  stochastic  methods  and  the 

deterministic  methods.  Stochastic  methods  are  those  methods  which  introduce 
random  models  of  uncertainty  and  assign  probability  measures  to  the  unknown 
parameters.  An  example  of  this  type  is  dual  control  [31].  In  this  approach, 
a system  is  simultaneously  identified  and  controlled  so  that  as  the  parameters 
are  estimated,  a control  is  implemented  in  such  a way,  that  for  these  estimates, 
the  system  behaves  in  some  ideal  sense.  As  the  estimates  are  made  to  converge 
in  a probabilistic  sense,  the  control  will  be  forced  to  converge  to  those 
values  that  elicit  the  ideal  behavior  in  the  system.  Deterministic  methods, 
on  the  other  hand,  use  mathematical  concepts  such  as  stability  and  algorithmic 
minimization  to  force  adaptation.  An  example  of  this  type  is  model  reference 
adaptive  control  (M.R.A.C.),  which  will  now  be  discussed. 
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CHAPTER  2 

THE  M.R.A.C.  METHOD  AND  THE  STABILITY  APPROACH 

In  the  Introduction,  it  was  suggested  that  many  control  problems 
subject  to  parameter  perturbations  may  be  solved  using  a M.R.A.C.  But  what 
is  a M.R.A.C.  and  how  is  it  realized?  To  answer  these  questions  we  now 
review  some  standard  definitions  and  results  which  may  be  useful  in  under- 
standing the  underlying  concepts  of  M.R.A.C.  and  also  make  clearer  the 
difficulties  of  the  methods  involved  [3] . 

The  M.R.A.C.  system  may  be  defined  as  an  adaptive  control  system 
which  uses  an  external  model  as  a reference  for  measuring  the  performance 
index  of  the  plant  (see  Figure  1).  It  will  be  shown  that  this  method  may  be 
employed  to  force  the  plant  responses  to  follow  the  responses  of  the  reference 
model. 

In  the  system  structure  of  Figure  1,  the  reference  model  may  be 
real  or  simulated.  The  error  vector,  representing  the  deviation  between  the 
real  and  ideal  responses,  is  defined  as  the  difference  between  the  plant  and 
model  outputs.  Generally,  this  vector  may  include  any  r all  of  the  state 
variables  along  with  other  outputs.  Once  the  "performance  measures"  are 
taken,  the  errors  are  then  fed  back  through  an  adaptive  algorithm  which  will 
either  generate  parameter  perturbations  or  alter  the  input  command  u, 
depending,  respectively,  upon  whether  the  parameters  are  physically  available 
for  change  or  not.  It  will  be  this  adaptive  loop  which  forces  the  error 
vector  and  its  derivative  to  approach  zero  in  the  case  of  perfect  model 


reference  adaptation. 


Figure  1.  M.R.A.C.  Structure 
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originating  with  Rucker  [14],  state-variable-filters  (SVF)  exploited  the 
commutation  of  linear  operators  in  linear  systems  and  thus  were  able  to 
successfully  measure  equivalent  state  variables  without  differentiation.  This 
was  indeed  a very  important  step  since  differentiation,  the  alternative  to 
SVF,  has  adverse  affects  on  noisy  measurements.  The  SVF  concept  was  then 
coupled  with  the  Kalman-Yakubovich  Lemma  [ 15]  by  Monopoli  [ 16] , Narendra  [ 10] , 
and  others  to  complete  the  solution  of  the  L.T.I.  M.R.A.C.  in  1972-3.  The 
additional  effects  of  the  Lemma  were  to  guarantee  a solution  for  the  L.T.I. 
problem  when  the  linear  reference  model  was  asymptotically  stable  and  also 
reduce  the  error  vector  stability  equation  to  a set  of  scalar  equations 
defining  the  adaptation  algorithm. 

The  hyperstability  approach  originated  in  1969  when  Landau  [17] 
first  proposed  it  as  an  alternative  for  the  solution  of  the  L.T.I,  M.R.A.C. 
problem.  Landau  later  generalized  the  class  of  potential  adaptive  laws  [18] 
and  then  with  the  aid  of  SVF  almost  completed  the  L.T.I.  case  design  [19]. 

A general  solution  for  generating  a positive  real  function  given  totally 
unknown  parameters  is  as  yet  unknown. 

The  M.R.A.C.  for  nonlinear  plants  represents  an  additional  level 
of  complexity,  hence  the.  solution  has  developed  more  slowly.  One  of  the 
earlier  contributions  was  published  by  Lindorff  and  Monopoli  [20]  in  1967. 
Using  Lagrange  stability  theorems,  they  produced  a solution  requiring  restric- 
tions on  the  input-output  relations,  knowledge  of  the  parameter  and  parameter 
derivative  bounds,  and  all  state  variable  measurement.  The  solution  was  then 
broadened  by  Lai  and  Mehrotra  [21]  in  1972  using  an  absolute  stability  theorem 
by  Davison  [22].  This  solution  required  state  variable  measurements  and 


Lj ^ 
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restricted  the  formulation  for  the  reference  model,  but  it  allowed  a very 
general  calss  of  nonlinearities  in  the  plant  and  no  parameter  bound  requirements. 
Lowe  and  Rowland  [23]  then  introduced  a Lyapunov  stable  M.R.A.C.  which  required 
parameter  and  parameter  derivative  bounds,  but  allowed  very  general  nonlinearities 
and  simplified  the  adaptive  control  to  the  bang-bang  type.  Finally,  the  most 
recent  major  improvement  came  from  Monopoli  and  Hsing  [24]  in  1975  when  they 
introduced  the  SVF  concept  into  the  nonlinear  problem.  This  was  a welcome 
result  since  it  reduced  the  need  for  differentiation.  However,  to  apply  the 
filters,  they  had  to  make  several  restrictions  as  discussed  below. 

It  was  shown  by  Monopoll  and  Hsing  [24]  that  SVF  may  be  used  to 
reduce  the  need  for  all-state-variable  measurement  in  a certain  class  of 
plants.  This  class  of  plants  is  of  the  following  form: 

i 

D.(p)x(t)=  S N (p)  u (t)  + c f . (x  , . . .X  t ) , l*'l,...,q 

1 i ij  J 1 1 i q 

where  Uj(t),  j=l  to  / and  x^,  i=l  to  q are  the  plant  inputs  and  outputs 
respectively;  f^(Xj^, . . . ,x^,  t)  are  nonlinear  time-varying  functions  of  known 


form;  p is  the  differential  operator  d/dt; 


m. 

1 


m 


i-1 


Di(P)-P  +.iiP 


J 


m 


N, .(p)  =b..p  +b..  p ... 


ij 


+b 


ijm 


i-1,, 


ij 


^i"  ^^ii  *^i  unknown  constant  or  slow-time-varying  parameters;  plus 

other  assumptions. 

Clearly,  under  the  above  conditions,  all  the  state  variables  are 
functions  involving  linear  operations  on  the  outputs,  and  the  are  not 
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functions  of  derivatives  in  x . Therefore  linear  operators  (i.e.,  SVF)  may 

q 

be  used  to  recover  the  states  from  the  outputs  without  derivatives.  However, 
when  the  above  assumptions  are  violated,  SVF  are  of  little  use  since  linearity 
can  no  longer  be  exploited.  Therefore,  due  to  lack  of  a better  technique, 
plant  states  must  generally  be  measured. 

The  above  development  shows  that  several  stability  synthesis  tech- 
niques have  been  used  to  treat  the  nonlinear  case.  The  Lyapunov  and  hyper- 
stability  methods,  however,  appeared  to  be  the  most  general.  Hence,  we  now 
review  some  standard  definitions  and  theorems  with  the  intention  of  outlining 
and  comparing  the  two  methods . 

m 

Definition  [5] : A function  V(x),  V:R~*R,  xeR  is  a Lyapunov  function  if  it 

satisfies  the  following  conditions: 

1)  V(x)  is  continuous  and  has  continuous  first  partial  derivatives 
in  some  open  region  R about  the  origin,  defined  by  l|x||  < a, 
where  a is  a positive  constant. 

2)  V(x)  is  positive  definite  in  the  region  R, 

3)  Relative  to  the  system  x“F(x),  F a function  vector,  the  time 
derivative  V(x)  along  the  trajectories  of  the  system  is  negative 
semlde finite  in  R. 

Theorem  (Lyapunov)  [5]:  If  for  the  system  x*F(x),  F(0)“0,  a Lyapunov  function 

can  be  found  in  some  neighborhood  R about  the  origin,  then  the  origin  is 
stable . 
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Theorem  (Lyapunov)  [5] ; For  the  system  x»F(x),  F(0)*-0.  If  in  a neighborhood  R 
about  the  origin  a Lyapunov  function  V(x)  can  be  found  such  that  V is  negative 
definite,  then  the  origin  is  asymptotically  stable. 

The  following  will  be  used  to  define  hyperstability  [6] : Consider 

the  system  described  by  the  state  equations  of  the  form 

X = F X + Gu 

y - H X + Ju  (1) 

where  the  triple  (F,G,H)  is  observable  and  controllable,  and  u and  y have 
the  same  dimension.  Define  u(*)  such  that  for  all  T 
T 

1 u' (t)y(t)dt  ^ f>[lix(0)ll]  sup  lix(t)ll  (2) 

o o < t < T 

where  6 is  a positive  constant  depending  on  the  initial  state  x(0)  of  the 
system  but  independent  of  the  time  T and  the  norm  is  the  Euclidean  norm. 

Definition  (Popov)  [6] : The  system  (1)  is  termed  hyperstable  if  for  any 

u(*)  in  the  subset  defined  by  (2),  the  following  inequality  holds  for  some 
positive  constant  K and  all  t: 


llx(t)ll  < k(llx(o)ll  + 6) 


(3) 


If  in  addition  u(*)  is  bounded  and  lim  x(t)  * 0,  the  system  fs  termed  asympto- 

00 

tically  hypcrstable. 


Definition  [ 7] ; A rational  fractional  function  M(s)  of  the  complex  variable  s, 
real  valued  for  real  values  s is  called  positive  real  if 


1)  Re  M(Jw)  >0  V W e R , 

2)  it  has  no  right-hand-plane  poles,  and 

3)  it  has  single  imaginary  poles  with  real  and  positive  residues. 
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Definition  [6] : A rational  transfer-function  matrix  Z(s)  is  termed  positive 

real  if  the  following  three  conditions  are  satisfied: 

1)  Z(s)  has  real  elements  for  real  s. 

2)  The  elements  of  Z(s)  are  positive  real. 

3)  For  any  real  value  of  W such  that  no  element  of  Z(jw)  has  a 

f 

pole  for  this  value,  z(jw)  + z (jw)  is  positive  semidefinite 
Hermit ian. 

Theorem  (Popov)  [6] : The  system  defined  by  (1)  and  (2)  is  hyperstable,  if 

and  only  if  the  transfer  function  of  (1)  is  positive  real. 

Theorem  (Popov)  [6] : The  system  defined  by  (1)  and  inputs  (2)  is  asymptotically 

hyperstable,  if  and  only  if  the  transfer  function  z(s)  = H(SI-F)'^G+J  for  the 
system  (1)  is  strictly  positive  real. 

Using  the  Lyapunov  theorems  and  the  Lyapunov  functions  together 
with  some  algebra,  it  was  shown  by  Parks  [8],  Monopoli  [9],  and  Narendra  [10] 
that  the  M.R.A.C.  problem  reduces  to  synthesizing  a Lyapunov  function  such 
that  if  the  error  vector  of  Fig.  1 is  chosen  as  the  state,  the  associated 
error  system  is  asymptotically  stable.  This  synthesis  generally  involves 
the  introduction  of  a Lyapunov  functi^'n  V which  is  quadratic  in  the  error 

vector,  forming  the  associated  V and  then  manipulating  parameter  perturbation 

• • 

signals  in  V such  that  V is  negative  definite.  Thus,  the  generation  of  the 
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V function  and  the  manipulation  of  the  parameter  signal  perturbations  in  V are 
the  key  to  the  solution  using  the  Lyapunov  approach,  both  of  which  are  often 
difficult  in  general. 

In  hyperstability  synthesis,  on  the  other  hand,  one  begins  by 
restructuring  the  system  of  Fig.  1 so  that  the  error  vector  is  defined  as  the 
state  and  a partionlng  of  forward  and  feedback  system  signals  yields  a linear 
forward  path.  One  then  designs  a constant  matrix  D such  that,  if  D operates 
on  the  forward  path  signals,  the  overall  forward  transfer  function  is  positive 
real  (see  Fig.  3)*,  The  synthesis  is  completed  by  finding  functionals  for 
the  forward  signals  such  that,  when  fed  back  through  the  feedback  path,  the 
input-output  relations  of  the  feedback  path  satisfy  an  inequality  analogous 
to  (2).  Hence,  the  M.R.A.C.  problem  is  reduced  to  a three  step  process,  the 
last  two  steps  being  difficult  to  do. 

Since  Lyapunov  and  hyperstability  synthesis  are  both  difficult  to 
apply  and  have  equally  desiable  stability  properties,  the  choice  depends  upon 
the  particular  application.  In  this  paper  the  hyperstability  was  chosen  since 
few  results  existed  for  the  hyperstable  nonlinear  M.R.A.C.  and  more 
importantly  because  it  appeared  to  be  the  simpler  design. 


The  new  structure  is  called  the  Equivalent  Popov  Structure. 
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TABLE  1 


SUMMARY  OF  LYAPUNOV  AND  HYPERSTABILITY 


DESIGN 


LYAPUNOV 


HYPERSTABILITY  METHOD 
USED  HERE 


Form  error  vector  between  model 
and  plant  states. 


1.  Same  as  Lyapunov 


2.  Define  the  error  vector  as  the 
state  of  the  system  and  algebrai- 
cally manipulate  plant  and  model 
state  equations  to  form  system 
state  equations  of  the  form 

& - f(e,x  ,x  ,u,t) 
n p 

3.  Find  a Lyapunov  function  and  form 

V. 


2.  Same  as  Lyapunov 


3.  Choose  M<  as  the  sum  of  variable 

gains  times  the  functions 

of  X and  x appearing  In  the 
m p 

plant  and  model. 

4.  Substitute  u from  step  3 Into 

f of  step  2 and  separate  f Into 
linear  and  nonlinear  parts. 


4.  Try  to  define  u as  a function  of 

X ,x  and  t such  that  V Is  negative 
m*  p ° 

definite.  If  such  an  u exists, 

design  Is  finished.  It  not,  go 

back  to  step  3. 


Define  a matrix  D such  that  D 
times  the  transfer  function  of 
the  linear  part  found  in 
step  4 Is  positive  real. 


Try  to  find  variable  gains 
of  step  3 as  a function  of 
x^,x  such  that  input -output 

relations  of  the  nonlinear  part 
found  in  step  4 satisfy  an 
Integral  inequality  analog  to 
(2).  If  this  step  and  step  3 are 
satisfied,  solution  is  complete. 
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CHAPTER  3 


SOLUTION  FOR  THE  NONLINEAR  M.R.A.C,  PROBLEM 


USING  HYPERSTABILITY 


In  this  chapter  It  will  be  shown  that  the  design  method  of 
multivariable  adaptive  model  following  control  described  by  Landau  [19] 
may  be  extended  to  Include  a class  of  nonlinear  multivariable  plants. 
ASSUMPTIONS ; 

1.  The  nonlinear  plants  to  be  treated  here  are  of  the  form; 

X - A X + E r“  f“(x  ) + B u (4) 

P P P a.l  P P P P P 

where  x Is  an  m-dlmenslonal  state  column  vector  available  from  noiseless 
P 

a 

measurements:  u Is  an  r-dlmenslonal  Input  column  vector;  A and  R for 
P P P 

or~l,...,3  are  mxm  unknown  constant  matrices;  B^  Is  an  n^r  constant  matrix 

whose  possibly  unknown  elements  b^^  which  must  lie  within  bounds  defined 

a 

by  the  real  numbers  p^^  and  for  o?-l,,..,P  Is  an  n-dlmenslonal 

column  vector  of  known  nonlinear  functions  of  the  state  x^;  the  plant 

(4)  satisfies  the  existence  and  uniqueness  properties  for  a solution 

(l.e..  Global  Llpschltz  conditions  [33]);  and  P ■ greatest  Integer  of  [k/m] , 

where  k Is  defined  as  the  minimum  number  of  functions  f^  needed  to  represent 

each  of  the  plant  nonllnearltles  F,  in  the  form  F.  ■=  .^,a..f.  with  a.  .'s 

j j 1«1  Ij  1 Ij 

elements  of  the  real  numbers. 


An  example  for  the  last  condition  Is: 
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^1 

“ ®11^1 

+ 8^3X3  + 

r(x2X^  + cos  Xj^)  + 

x„ 

■ 

+ r„„  cos 

X- 

2 

22  2 

22 

1 

X3 

“ ®33’^3 

+ ’’32"2 

Here  P ■ greatest  Integer  of  [2/3],  implies  P*1  and 


®11  ° ®13 

r r 0 

X2X3 

’^l  ° 1 

1 

0 a22  0 

X + 

P 

0 r^2  0 

cos  Xj^ 

+ 

{ 

0 0 

1 

0 0 a^^ 

0 0 0 

0 

• — 

^2_ 

Uii 


Essentially  we  are  considering  nonlinearities  f which  factor  with  respect 
to  a set  of  parameters.  That  is,  for  parameter  vector  P,  f(x,P)  ■ PF(x) 
where  F is  a column  vector  of  appropriate  dimension. 

2.  The  model  structure  will  be  given  by 


X * A X + 2 F^(x  ) + B ,u 

m mm  Y“1  m m m mm 


(5) 


where  x is  an  m-dimensional  state  column  vector;  u is  an  A -dimensional 
m m 


command  input  column  vector;  is  a known  nXn  Hurwitz  matrix;  and 


for  \ = 1,...,P  are  nX^  and  nXn  known  constant  matrices  respectively; 


Y 

F (x  ) for  Y=1,...,P  is  an  n-dimenslonal  column  vector  of  the  nonlinear 
m m' 


functions  of  the  states;  and  P is  defined  similarly  to  P but  with  respect 


to  the  model  nonlinearities. 
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3.  The  plant  and  model  satisfy  perfect  model  following  conditions. 

"Model  following  Is  a method  of  control  In  which  the  plant  Is 

controlled  to  behave  like  an  Ideal  model."  When  It  Is  "possible  to  match 

perfectly  the  dynamics  of  the  compensated  plant  with  those  of  the  model" 

It  Is  called  perfect  model  following  or  matching  [25].  Hence,  perfect  model 
following  conditions  are  those  conditions  Insuring  the  existence  of  a feedback 
compensating  control  such  that  the  plant  and  model  are  Input-output  equivalent. 

Under  some  circumstances  (Chan  [25]  for  the  linear  case).  It  can 
be  shown  that  when  perfect  model  following  conditions  are  not  satisfied,  the 
error  may  still  be  bounded.  Therefore  without  prior  knowledge  of  the  plant, 
perfect  model  following  may  still  be  useful.  Sufficient  conditions  for  perfect 
model  following  for  the  class  of  nonlinear  plants  treated  here  will  be  derived 
In  Appendix  B. 

Comment : Due  to  the  distributive  nature  of  the  nonllnearltles  and  the 

technique  to  be  snown,  It  will  only  be  necessary  to  show  the  solution  for 
3“1  and  p«l.  Hence,  It  will  be  assumed  that  3=1  and  P=1  In  all  that  follows 
unless  stated  otherwise.  The  treatment  for  3>1  and/or  P>1  will  become  obvious 
as  we  proceed. 

STRUCTURE  OF  THE  CONTROLLER: 

There  are  two  basic  adaptive  model  following  structures,  signal 
synthesis  and  parameter  adaptation.  While  the  two  are  mathematically  equiva- 
lent [19],  signal  synthesis  provides  adaptation  by  modifying  the  feedback 
gains  Instead  of  directly  changing  the  parameters  as  Is  the  case  In  the  latter. 


since  usually  not  all  the  parameters  are  available  for  adaptation,  signal 
synthesis  is  the  natural  choice.  Figure  2 shows  how  the  signal  synthesis 
structure  is  used  in  our  particular  case. 
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I 

I 

I 

I 

I 

I 
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Having  chosen  signal  synthesis,  we  define 


A ^ 2 

u = u + u 
p p p 


1 A 

u =-Kx  +Kx  +K„F  +K„„F  +KR 
p p p mm  ^ 


(6) 


u ^ = AK  X + AK  F + AK  F + AK  u 
p pp  Rp  Rm  urn 


v;i»crc  is  the  plant  control  signal;  u^  is  some  nominal  control  whose  design 

2 

shall  be  discussed  later;  u'  is  the  adaptive  control;  K , K , K , K , AK  , 

p p m m Rp  p 

AK  , and  are  rXn  weighting  matrices;  and  and  AK^  are  also  rXX  weighting 
m 

matrices , 

There  are  two  reasons  for  choosing  the  feedback  signal  input 
structure  (6).  First,  recall  that  the  parameters  factor  out  of  the  known 
nonlinearities.  Hence,  it  appears  that  under  certain  conditions,  parameter 
adaptation  may  be  used  to  add  and  subtract  functions  in  the  plant  state 
equation  by  feeding  back  the  required  functions  and  states.  Justification 
for  the  problem  treated  appears  later  and  in  Appendix  li . 

The  second  reason  for  choosing  the  structure  (6)  is  that  it  may  be 
desirable  to  bias  the  input  control  level.  For  example,  if  circuit  implemen- 
tation places  magnitude  constraints  upon  the  adaptive  gains  and  rough  estimates 
are  available  for  the  parameters,  conventional  feedback  may  be  used  to  drive 


I 


J 


f 
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the  system  within  the  range  of  the  adaptive  gains.  In  this  problem,  one 

biases  the  control  by  defining  as  a coarse  nonadaptive  control  stabilizing 

2 

the  system  for  the  assumed  nominal  values  of  the  plant.  Likewise  u^  is 
defined  as  a fine  adjustment  adapting  to  any  perturbations  from  the  assumed 
values . 


It  should  be  stressed  that  the  assumed  nominal  values  need  not  be 
close  to  the  actual  values,  but  the  closer  they  are,  the  better  the  biasing 
effect . 


PERFECT  MATCHING  AND  THE  EQUIVALENT  POPOV  STRUCTURE: 

Having  chosen  our  structure  and  controls,  one  must  ask  under  what 
conditions  will  the  control  force  the  plant  to  follow  the  model.  This 
question  is  partially  answered  in  Appendix  B,  where  it  is  shown  that  sufficient 
conditions  for  perfect  model  following  are  that  the  error  vector  e equal 
zero  and  that  equation  (Bl)  be  satisfied.  It  is  also  shown  that  control 
law  (B2)  will  work  under  the  above  conditions. 

Prior  to  using  the  Assertion  of  Appendix  B to  formulate  the  nominal 

control  u^  , we  must  first  assume  that  (Bl)  is  satisfied.  Justifying  this 

point  is  difficult  since  it  is  impossible  to  tell  a priori  whether  B^  is 

non-singular  or  whether  the  matrices  (A  -A  ),  B , R and  R lie  in  the  null 

m p m m p 

space  of  (I-BpBp).  However,  given  some  prior  information  we  may  be  able  to 
deduce  that  (Bl)  is  true.  If  not,  we  may  take  comfort  in  remembering  that 
(Bl)  is  only  a sufficient  condition.  Also  Chan  [25]  has  shown  that  even 
without  perfect  following  our  error  may  still  be  bounded  and  thus  the  control 
may  still  be  acceptable. 
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Accepting  (Bl)  to  be  true,  we  choose  K , K , K , K , and  IC 

m p u Rp  Rjn 

by  setting  the  A K's  = 0 and  choosing  nominal  values  for  the  parameters  in 
(B2).  That  is, 


K-K  = B(A-A);K  = B B; 
m p p'mp  v pm 


K = -B  R 
Rp  P P 


; K_  = B R ; 
K-m  P 


where  the  means  nominal. 

Note  that  the  AK's  will  be  designed  to  compensate  for  any  errors 
in  the  assumed  nominals  and  setting  K^“0  does  not  affect  the  controller. 

Now  consider  the  design  of  the  A K's.  We  will  construct  an 
equivalent  Popov  structure  and  then  choose  adaptive  laws  in  such  a way  that 
will  make  the  overall  system  asymptotically  hyperstable. 

The  selection  of  the  Ak's  will  be  accomplished  using  the  hyper- 
stability concepts  that  were  described  earlier.  The  M.R.A.C.  structure  of 
Figure  2 will  be  restructured  to  an  equivalent  Popov  structure.  This 
restructuring  is  not  difficult  to  do  and  forces  the  M.R.A.C,  structure  of 
Figure  2 to  be  characuerized  by  a system  having  a linear  forward  path  and 
nonlinear  feedback  path  in  the  error  system. 

The  following  procedure  will  produce  an  equivalent  Popov  structure 
for  system  described  by  (4),  (5),  and  (6) : 

Define  the  error  as 


Subtracting 


A 

e * X - X 


m e 

(4)  from  (5)  and  substituting  (6)  and  (7)  yields 


(7) 
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e=Ax  +RF  +Bu  - Ax  -RF  - B[(AK-K)x 
mm  mm  mm  pp  pp  p ppp 

Vo,  * %,  ) ‘'m  '‘'u 

p p mm 

Now  adding  and  subtracting  terms  on  the  right 


e = (A  - B K ) e + B (b"^  R - K„  - AK  ) F + B [ b'''b  - K - AK  ] 

^ m pm'  p^pmR  R^m  PPm  u u 

mm 


u + B [b‘*’(A  - A)+K  -AK  - K]x 

~ pp  m p'  p p mp 


m 


(8) 


-B[B  R +IC  +AK„]F 
P P P \ \ V 

Then  separate  (8)  into  linear  and  nonlinear  parts  by  defining 


" ■ -»1  ■ '“St  + % - “p  “J''. 

m m 


+ (K  + AK  - B B )u 
u u p m m 


+ (K  + Ak  - K - B (A  - A )]x 
m p p p m P P 


+ |BpEp  + (AK^  +K^))Fp 


(9) 


and  substituting  (9)  into  (8)  we  obtain 


e = (A  - B K )e  + B W, 
m pm  pi 


(10) 


Define  the  following  linear  operator 


V = De 


(11) 


Therefore  equations  (9),  (10)  and  (11)  establish  the  equivalent  Popov  structure 
for  the  system  (4),  (5)  and  (6),  See  Figure  3. 
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Using  the  definition  of  hyperstability  given  earlier,  we  reiterate 
that  the  system  described  by  (7),  (9),  (10),  and  (11)  is  hyperstable  if  there 
exists  constants  6 > 0 and  y > 0 such  that  any  solution  of  the  system  satisfies 
the  inequality 

lle(t)ll  < 6(lle(C)ll)  + y ; Vt>0;  6>0;  y>0  (12) 


for  all  w and  v such  that 


Tl(0,t)  “ f w'^V  dT  > -y^  - y^  ; 


sup 

o < y ^ t 


|e(T)||y^  > 0 


(13) 


and  is  asymptotically  hyperstable  if  in  addition  to  the  above 
for  all  W and  V satisfying  (13). 


lim  e(t)«0 
^^00 


Therefore,  by  the  Popov  theorem  necessary  and  sufficient  conditions 
that  the  system  described  by  (7),  (9),  (10),  and  (11)  be  asymptotically 
hyperstable  is  that  the  transfer  function  matrix  characterizing  the  linear 
part  of  the  equivalent  Popov  structure, 

Z(S)  - DiSI  -(A^-BpK^)]  , (14) 

be  strictly  positive  real. 

Thus  we  can  force  e-*o  by  searching  for  a constant  matrix  D to 
make  Z(s)  strictly  positive  real  and  by  choosing  the  Ak's  in  (9)  in  such  a 
way  that  (13)  is  satisfied. 

THE  ADAPTATION  RULES  AND  THE  MATRIX  D: 

The  adaptation  rules  are  the  adaptation  algorithms  that  will  guarantee 
the  error  system  is  asymptotically  hyperstable,  and  the  following  assertion 
provides  such  an  algorithm. 
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Assertion:  The  following  choices  for  the  Ak's  satisfy  (13): 


AKj^  - J LV[SFj^  dT  + LV(SFj^ 


AKj^  - J I)y[QF  dT  + 0V[QF j'^  + AK^^  (0) 
p o P P p 


^ m fn 

Ak  “ J NV[PX  dT  + NV[PX  ]^  + AK  (0) 


Ak  = f MV[Tu  dT  + MV[Tu  ] ^ + AK  (0) 

II  «J  m m II  ' ' 


where  L,L,S,M,M,T,N,N,P,0,0,  and  Q are  constant  positive  definite  symmetric 


matrices . 


This  assertion  is  the  most  important  part  of  this  approach  and 


is  a direct  extension  of  the  Landau  proportional+integral  control  to  include 
the  system  given  by  (4).  A proof  of  the  assertion  is  given  in  Appendix  A. 

The  object  of  the  matrix  D is  to  make  Z(s)  in  (14)  strictly  positive 
real.  However,  as  can  be  seen  from  (14),  some  upper  and  lower  bounds  for 
each  of  the  parameters  of  are  needed  to  insure  the  positivity  requirement. 

Under  these  conditions  tests  for  positive  realness  such  as  those  of  Anderson 
[26]  or  Silyak  [27]  may  prove  to  be  useful. 

In  the  special  case  when  Bp  is  known  completely,  but  other  parameters 
in  the  plant  are  unknown,  the  following  theorem  of  Kalman  and  Popov  provide 
a satisfactory  solution  for  D.  See  Landau  [19]. 

Theorem  (Kalman,  Popov)  [19]:  The  Hurwitz  matrix  A and  matrices  b,C  and  J 

define  an  observaljJe  and  controllable  . system; 

A theorem  derived  by  landau  [ 18]  may  be  used  to  generalize  these  matrices 
to  include  positive  real  kernels.  See  Appendix  A. 
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k\  - J + 
ta  + A^r  = 

B^r  + - C (16) 

In  our  case  Z(s)  ■ D(SI  -A  +BK)^B  with  J*0.  Therefore, 

m pm  p 

if  (A,  - B K ) is  Hurwitz,  then  B^  T ■ D will  satisfy  the  above  theorem  for 

m p m p 

positive  realness.  Note  that  if  (A  - B K ) is  not  known  to  be  Hurwitz,  K 

m p m m 

may  be  set  to  zero  as  mentioned  earlier  and  thus  T will  exist.  (Recall  A^ 
was  postulated  to  be  Hurwitz  in  the  problem  definition.) 

At  this  point  the  design  procedure  has  been  completed  and  all  that 
remains  is  to  illustrate  the  theory  by  using  it  on  a particular  problem. 

Ir.  Chapter  4 we  apply  the  theory  to  a specific  nonlinear  plant  and  simulate 
the  resulting  system. 


I 

I 

\ 

) 

1 

I 
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CHAPTER  4 
AN  EXAMPLE 


In  the  last  chapter,  a design 
problem  was  given  in  general  terms.  In 
to  a specific  nonlinear  plant  that  will 
model . 

The  Problem; 

The  plant 


w 


w 


w 


1 

2 

3 


^2~^3 
I^  ”2”3 

^3‘^1 


^l'^2 
I.  ’'l''2 


procedure  for  the  nonlinear  M.R.A.C. 
this  chapter  we  will  apply  the  procedure 
be  forced  to  adapt  to  a linear  uncoupled 


u 

1 


models  the  tumbling  motion  of  an  orbiting  satellite  whose  tumbling  rate  far 
exceeds  the  orbiting  rate  (Hsu  and  Meyer  [28])  • The  W's  are  angular  velocities, 
I's  moments  of  inertia,  and  the  U's  are  torques  due  to  thrust  rockets.  Note 
that  I^  ^ 0 (i  = 1,2,3).  Therefore  is  nonsingular,  and  thus  the  perfect 
following  conditions  (Bl)  are  satisfied. 

Assume  for  the  purposes  of  the  example  that  the  actual  values  of  the 
parameters  1^,  I2,  and  I^  are  10,  40,  and  25  respectively.  Thus  the  plant 
is  given  by 


”1 

1.5  0 o' 

"2  "3 

''2 

s 

0 375  0 

”1  ”3 

+ 

D 

0 0 -1.25 

Wi  w^ 

D 

-■P 


(18) 
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Also  assume  the  only  information  given  about  the  plant  is  its 

structure  (i.e.,  none  of  the  parameters  are  known  exactly)  and  the 

following  bounds  for  the  R and  B matrices  which  may  have  been  found  through 

P P 

experimentation.  (For  this  plant,  these  bounds  represent  a 10:1  range  over  the 
actual  values). 


Rp~  (-15,15) 


(-15.15) 


(-15,15) 


(0,lj  0 0 

0 (0,1]  0 

0 0 (0,1] 


Given  (19),  we  may  choose  nominal  values*  for  R^  and  anywhere 
in  the  given  range.  For  convenience  choose 


R„  “ [0]  and  - [l] 

P P 


Suppose  the  reference  model  is: 


+ 0 


-1  w., 


1 u 


THE  DESIGN: 


We  now  set  the  Ak's  in  (B2)*’'  to  zero  and  substitute  the  nominal 


values : 


As  stated  in  Chapter  3,  choice  of  nominal  values  only  affects  the 
separation  of  magnitude  scales. 


Since  it  has  been  shown  that  K may  be  set  to  zero,  we  will  do  so  for 
convenience. 


Since  A does  not  appear  in  the  plant,  it  is  set  to  zero.  Hence 


The  operator  D is  now  chosen  to  make  Z(s)  in  (14)  positive  real 
We  note  that  for  our  case 


S+1 

0 

0 

-1 

’ (0,1] 

0 

0 

0 

s+1 

0 

0 

(0,1] 

0 

0 

0 

s+1 

0 

0 

(0,1] 

"(OJ-J-  0 0 1 

Z(S)  - D 
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and  therefore  choosing  D ■ [l]  assures  that  Z(S)  Is  strictly  positive  real. 

Assign  for  simplicity  the  adaptive  gains  L“L*M"M*N“N" 
Q * 0 equal  to 


The  Ak's(O)  ■ [0] , and  set  S ■ T ■ P ■ q equal  to  the  value  in  (21) 


" «im  - "ip  * 3 


V “ e 
i 1 


u ■ u + u 
P Pi  P2 


u ■ - w.  + w. „ 
p,  . ip  im 


u ■ [ r 100  v,w,  dT  + 100  v.w,  ]w, 

P2i  i 1 IP  ^ 1 IP  IP 


+ [J  100  v^W2pW3pdT  + 100  ViW2pW3p]w2pW3p 

+ [f  100  v,u  dT  + 100  v,u  ]u 
•J  1 m,  1 m,  m, 

o 1 11 


Up  - [|  100  v^w^pdT  + 100  v^w^plw^p 

22  j 


+ [/  100  V2W^pW3pdT  + 100  V2WjpW3p]w^pW3p 

+ [ f 100  v-u  dT  + 100  v.,u  ] u 

i 2 m2  2 m2  m2 


f 

f 

r 


I 

1 

I 

I 

I 

I 

1 

I 
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P23 


u 

■ [J  100  v^w^pdT  + 100  v^w^^Jw 


3"3p'"3p 
t 

+ [f  100  v_w.  w.  dT  + 100  v_w_  w ]w„  w 
o 3 2p  Ip  3 2p  Ip  2p  Ip 

+ [J  100  v-u_dT  + 100  v-u_]u 

_ j J j j m_ 


(26) 


We  now  present  the  results  and  problems  encountered  with  the 
simulation  of  the  adaptive  system  described  above. 

SIMOIATION  AND  RESULTS 

The  main  objective  of  the  simulation  is  to  show  that  adaptation 
worked  with  and  without  initial  conditions  in  the  plant.  We  were  to  confirm 
this  by  measuring  the  error  vector  between  model  and  plant  with  and  without 
the  adaptive  controller.  If  given  the  same  initial  conditions,  the  error 
vector  of  the  adaptive  system  converged  near  zero  and  non-adaptive  system 
didn't,  then  the  adaptation  was  working. 

Other  objectives  were  to  study  the  effects  of  the  adaptive  gain 
constants  on  error  transients  and  convergence  rates,  and  to  study  the  use 
of  plant  and/or  model  states  in  the  formation  of  the  nonlinearities  in  the 
adaptive  loop.  The  adaptive  gain  study  is  important  since  theoretically 
the  gains  may  be  arbitrarily  high  and  thus  produce  arbitrarily  fast 
adaptation.  The  other  study  was  introduced  since  it  can  be  shown,  by  an 

"if 

The  simulation  was  implemented  using  the  C.S.M.P.  language  on  a DEC-10 
digital  computer.  For  further  information  and  sample  simulation  programs 
see  Appendix  C. 

** 


Arbitrary  up  to  the  time  optimal  value. 
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argument  similar  to  that  used  in  chapter  4,  that  an  w.  w.  integral  may  be 

im  jm 

•)*r 

used  in  place  of,  or  in  addition  co,  the  Integral  in  equation  25 

and  still  have  an  asymptotically  hyperstable  error  system. 

The  non-adaptlve  plant  and  reference  mode  described  (equations 
18  and  21)  were  both  simultaneously  simulated  using  the  same  sinusoidal  Input 
and  zero  initial  conditions.  An  error  vector  was  formed  as  the  difference 
between  the  plant  and  reference  model  states,  and  recorded  as  f notions  of 
time  into  what  now  are  figures  4 to  6.  These  figures,  along  with  a 
typical  state  signal  of  the  plant  (Figure  7),  show  that  the  errors  are  very 
large  but  bounded. 

The  adaptive  system  (equations  18,  21,  and  26)  was  then  similarly 
simulated  with  the  same  inputs,  same  initial  conditions,  and  adaptive  gains 
of  10  (see  equation  25).  The  results  of  this  simulation  are  given  in 
Figures  8 and  9.  Comparing  Figure  4 to  Figure  8 and  Figure  6 to  9,  it  is 
plain  to  see  that  the  error  vector  for  the  adaptive  system  seems  to  be 
considerably  less  than  its  non-adaptive  counterpart  for  all  time,  and  that 
it  is  converging  to  within  a small  neighborhood  of  zero. 

Having  shown  that  the  system  adapts  for  zero  initial  conditions, 

the  system  was  then  simulated  with  various  initial  conditions.  Figures  10, 

11  and  12  depict  the  state  errors  converging  to  within  some  neighborhood 

of  zero  for  initirl  conditions  of  w°  = -3,  w°  •=  1 and  w”  - 4.  Initial 

Ip  2p  3p 


Integral  has  the  form 
t 

[ 100  V w.  w.  + 100  v V.  w.  ]w.  w. 

— im  jm  - im  jm  im  jt 


and  is  added  on  to  equation  25. 


32 

conditions  in  the  same  order  of  magnitude  produced  similar  results,  however 
when  the  initial  conditions  were  made  an  order  of  magnitude  higher  it  yielded 
inconclusive  results  (more  on  this  last  point  later  in  this  chapter). 

Figures  13,  14,  and  15  are  the  same  error  plots  of  the  Figures  10,  11, 
and  12  with  the  first  second  deleted  and  the  same  scale  blown  up.  Likewise  for 
Figures  16,  17,  and  18,  except  the  adaptive  gains  were  increased  to  lO/J  ~ 17. 
Comparison  among  these  figures  show  that  the  higher  gains  produced  faster 
initial  decreases  in  the  state  errors  by  a factor  of  at  least  1.5  to  2. 


However,  the  rates  leveled  off  after  about  10  seconds  to  produce  little 
improvement,  specially  for  w^p-  This  type  of  effect  has  been  observed  earlier 
while  simulating  an  adaptive  observer.  It  appears  that  high  gains  tend  to 


Increase  the  initial  improvements  but  induce  damped  tail  oscillations  as  the 

error  approaches  zero.  Quite  often,  decreasing  the  adaptive  gains  in  the  linear  i 


I 

I 

I 

I 

I 


terms  will  Increase  the  damping  rate  at  the  tall  end  at  the  expense  of  Initial 
rate  decreases. 

The  effects  of  replacing  by,  or  the  adding  of  an  term  were 

tested  by  similar  simulations  as  above.  Given  the  same  sinusoidal  inputs, 
adaptive  gains  of  10,  and  Initial  conditions  (w°p»3,  w°p  = -1,  and  w°p  = 0); 
systems  using  the  functional,  using  the  w^pWj^  functional,  and  using  both 

together  were  simulated.  The  results  are  shown  in  Figures  19  to  23.  Comparing 
Figures  19,  21  with  23,  and  20  with  22  we  see,  at  least  for  this  system,  that 
no  distinct  advantages  exist  for  any  of  the  cases.  That  is,  no  convergence 
properties  establish  one  case  better  than  the  other,  but  this  does  not  mean 


that  a noisy  version  of  the  system  will  exhibit  the  same  properties. 


* 

This  effect  was  noticed  while  simulating  an  adaptive  observer  using 
the  Narendra  structure  [30]. 
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Figure  8.  w^  plant  state  and  corresponding  performance  error  for  adaptive  system. 


Plant  Error 


Time  (Sec.) 


Figure  9.  W2  plant  state  and  corresponding  performance  error  for  adaptive  system. 
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System  signals  with  adaptation  and  Initial  conditions 
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Figure  10.  performance  error. 
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Figure  11.  v»2  performance  error. 
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Figure  12.  plant  state  output  and  corresponding  performance  error. 
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Printing  Interval  “ .2  sec. 
Inputs:  “ml  “ " “m3  “ 

Adaptive  Gains  ■ 10 


System  signals  with  adaptation,  initial  conditions,  and  adaptive  gains  equal 
10»  First  second  of  response  is  deleted  to  provide  a more  detailed  view  of 
transient. 


Time  (Sec.) 


Figure  13.  Scaled  Wj^  performance  error, 
trror  W, 
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Figure  14.  Scaled  Wj  performance  error. 
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Figure  15.  Scaled  w^  plant  state  output  and  corresponding  performance  error. 

Initial  Conditions:  w°p”3,  w°p--l  Printing  Interval -.  15  sec . 

inputs:  V“u^-u^-sin.5t 

Integration  Interval  “ .002  sec.  Adaptive  Gains  » 17 


inputs:  V“u^-u^-sin.5t 


System  signals  with  adaptation.  Initial  conditions,  and  adaptive  gains 
equal  17.  First  second  of  response  Is  deleted  for  more  detailed  view 
of  transient.  Compare  with  Figures  13,  14,  and  IS. 


Error  W, 


Time  (Sec.) 


2 3 3 e 7 3 1*  • 1*  9 19  • 9 IS  • 


Figure  16.  Scaled  w^  performance  error  with  Increased  adaptive  gains. 
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Figure  17.  Scaled  W2  performance  error  with  Increased  adaptive  gains. 
Wj  Plant  Error 


Time  (Sec.) 


Figure  18.  Scaled  w^  plant  state  and  performance  error  with  Increased  adaptive  gains 

. . , ^ o o _ , Printing  Interval  - .2  sec. 

Initial  Conditions:  w.  “3,  w„„  - -1  ^ ° 
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Figure  19.  w,  and  w.  performance  errors  with  w 's  functional  adaptation  only. 
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Figure  20.  Wj  plant  state  and  corresponding  performance  error  with  w *s  functional 

p 

adaptation  only. 


Initial  Conditions:  wf  "3,  w?  • -1  Printing  Interval  ■ .2  sec. 

ip  2p 


’'2p-° 


Integration  Interval  ■ .02  sec. 


inputs:  u„i-u^-u^-,ln.5t 

Adaptive  Gains  “ 10 


23.  System  signals  with  w 's  and  w 's  functional  adaptation. 
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Problems  Encountered 

Several  problems  appeared  during  the  simulation  phase  that  hindered 

the  evaluation  of  several  of  the  performance  checks.  It  appears  that  the 

C.S.M.F.  language  we  used  was  implemented  with  a single  precision  second 

order  Runge  Kutta  integration  routine,  and  it  was  introducing  truncation  and 

convergence  problems.  For  example,  while  simulating  for  Initial  conditions  of 

w°  “40,  w°  *-30  and  w^  “15,  the  routine  overflowed  with  an  integration 

Pi  P2  P3 

Interval  of  .000025  sec.,  underflowed  for  a .000020  sec.  value,  and  converged 
for  a .000021  value.  Clearly  Inconclusive  results  since  it  is  impossible  to 
tell  whether  the  system  adapted  or  whether  the  Integration  routine  truncated 
system  values  to  force  it  to  converge.  A related  problem  was  the  number  of 
iterations  needed  to  carry  out  the  simulation  assuming  everything  worked,  the 
integration  interval  had  to  be  .005  seconds  for  many  cases.  This  Implies  that 
a 20  second  run  would  need  4000  iterations  or  about  1%  minutes  of  CPU  time, 
clearly  very  expensive  for  a single  simulation.  Similar  problems  to  all  the 
above  were  also  encountered  while  increasing  the  adaptive  gains. 
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CHAPTER  5 

NEW  DIRECTIONS  AND  CONCLUSION 

PROBIZMS  UNSOLVED; 

To  make  the  design  procedure  in  Chapter  4 competitive,  several 
problems  must  first  be  solved.  In  particular,  a systematic  method  of  generating 
the  linear  matrix  D would  not  only  simplify  calculations  but  It  would  extend 
the  method  to  Include  unknown  matrix  parameters.  A clue  to  such  a result 
probably  lies  In  works  dealing  with  the  positive  real  functions  criterion 
[26,27].  Another  very  important  problem  is  the  generation  of  state  variable 
measurements  without  using  derivative  operations.  Though  such  a method  has  been 
found  for  linear  systems,  the  nonlinear  case  Is  still  open. 

A problem  of  practical  nature  is  that  of  formulating  a systematic  way 
of  choosing  the  adaptive  gains  of  the  controller  to  meet  performance  criteria. 
Such  a result  would  be  very  useful  in  designing  controllers  for  systems  with 
moderately  fast  time-varying  parameters  or  with  infrequent  discontinuous 
parameter  perturbations. 

Related  problems  of  an  algebraic  nature  are  the  formulation  of 
necessary  and  sufficient  conditions  for  perfect  model  following  in  the  nonlinear 
case,  and  the  existence  of  a minimal  parameter  set  structure  to  uniquely  define 
nonlinear  systems.  These  solutions  would  generate  a minimal  number  of  adaptive 
loops  and  guarantee  perfect  adaptation. 


With  respect  to  adaptation  time  transients. 

As  it  stands,  there  is  an  adaptive  loop  for  every  parameter  in  the 
nonlinearity  which  must  be  adapted  plus  minimal  number  for  the  A and  B 
matrices.  Thus,  possibly  too  many  to  be  practical.  ^ 
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l£  the  minimal  parameter  problem  described  above  is  solved,  the  following 

ij 

approach  may  prove  useful  in  the  structural  adaptation  problem  when  the  order 

1 

of  the  system  is  known.  i 

i 

Consider  the  class  of  nonlinear  functions  having  factorable  unknown 
parameters  and  is  also  analytic  in  the  region  of  question.  Then,  if  the  F 

^Qf  I 

matrix  of  (4)  is  assumed  to  be  composed  of  the  first  few  terms  of  the  nonconstant 

part  of  the  Taylor  expansion  terms,  the  AK  may  be  forced  to  adapt  (by  the  method 

p 

shown  in  Chapter  4)  to  the  coefficients  of  the  Taylor  expansions  of  the 
unknown  functions.  Since  the  functions  were  assumed  to  be  analytic,  it  would 
seem  that  the  model  following  error  should  at  least  be  bounded  by  a function  * 

involving  the  order  of  the  error  in  the  Taylor  expansion.  Thus,  the  possibility 
of  structural  adaptation  may  exist  by  a sli^t  extension  of  the  work  in  Chapter  4 
and  a minimal  parameter  structure. 

CONCLUSION: 

Parameter  adaptive  control  has  been  studied  as  alternate  approach  when 
conventional  controls  fail.  In  particular,  the  M.R.A.C.  structure  was  discussed 
in  conjunction  with  stability  synthesis.  Lyapunov  and  hyperstability  were  then 
presented  as  general  but  mathematically  complex  synthesis  methods.  Choosing 
hyperstability,  the  theory  was  successfully  applied.  Results  show  that  the 
approaches  of  Landau  [19]  and  Erzberger  [29]  may  be  combined  to  yield  a new 
hyperstable  M.R.A.C.  for  a class  of  nonlinear  systems. 

In  solving  the  problem,  two  useful  assertions  were  Introduced,  one 
providing  sufficient  conditions  for  perfect  model  following  and  the  other  an 
adaptive  algorithm  for  the  nonlinear  plant.  A design  example  was  then 
formulated  and  the  resulting  system  simulated.  The  results  confirmed  the 
effectiveness  of  the  approach. 
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APPENDIX  A 

DERIVATION  OF  ADAPTIVE  LAWS 


I 


The  following  assertion  provides  an  algorithm  for  generating  the 
adaptive  ^K's. 

Assertion:  The  following  choices  for  the  AK's  satisfy  (13). 

m o 


AKj^  - J Lv[sf^]'^dT  + LvCsf  + AKj^  (0) 


AKj^  •»  J bvCqf  ]^dT  + ovLqfp]^  + ak^  (O) 


P o 


AK  - r NvCpx  ]^dT  + NvCpx^]'^  + AK  (0) 

n « n o n ' ' 


P “o  - - p 


AK  = j MvL'To^l^dT  +MvrTP^]'^  + AK  (0) 


u " ■"  m 
o 


where  L,  L,  S,  M,  M,  T,  N,  N,  P,  0,  0,  and  q are  constant  positive  definite 
symmetric  matrices. 

Proof:  See  also  Appendix  II  of  f 19] . 

From  (9)  and  (13) 

^ T + 

iKo.t,)  . ,!  » Ik,,  -kak,,  -b  kJf„  bt 

o m m 


2 

+ f n'^Ck  + AK  - b"^  15  lu  d'r 
'J  u u p m m 

o 

+ J n'^[k„  + AK„  - dT 

m ppp'mp  p 


o 

t. 


(Al) 


I 


46 


Substituting  (15)  into  (Ic) 

t 


1 


T](0.t^)  = J ^N[Sf^3^  dT’  dT 


o o 
t. 


+ J n'^[J  M'ru  dT ' - B }u  dT 

^ ^ m o m 

o o 


+ NvCPx  dT*  -C  }x  dT 

« p op 


T ,„i 


o o 
t. 


+ OvCqf dT*  -D  }p  dT 

o o " ^ 


where 


+ J (N^LV)Cf^  S FjdT 

o 

+ J (N'^MV)[u^Tu^JdT 

o 

+ J^(N^NV)[x][  Px  ]dT 


o 

t. 


P P 


+ J (n^0V)[fJqf JdT. 
o PH 


= B (0)  ■ constant  matrix 

° ^ m m 

B “ B^  B - K - Ak  (0)  * constant  matrix 
o p m u u'  ' 

C * b"*^(A  -a  ) + K - K - AK  (0)  = constant  matrix 
o p'mp  p m p' 


(A2) 


Bp  Rp  - Kj^  ' matrix. 


In  the  case  where 


0>1  and  P > 1 we  would  have  similar  terms  to  these. 
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2 

Sufficient  conditions  that  Tl(0,tj^)  ^ -Y^  - Yq  are  that  each  term  of  (2c)  be 
2 

greater  than  where  i ranges  over  the  number  of  terms.  Thus  the  last 

four  terms  of  (2c)  satisfy  this  condition  since  L,  S,  M,  T,  N.  P,  0,  and  Q 
are  all  positive  definite.  Hence,  it  remains  to  be  shown  that  the  terms  of 


1 


and  therefore 


c ’ 19 

J Y(T)[J  Y(T')dT']dT  - Y(T)dT]^  . 


o o 


Hence  (4c)  equals 


which  is  obviously 


ifju;  -*oiP 


It  is  now  possible  to  show  that  the  positive  definite  matrices  in  (15)  may  be 
generalized  to  the  class  of  positive  definite  kernels  [18]. 

Definition:  The  square  matrix  kernel  K(t)  is  termed  "positive  definite"  if 

for  each  interval  [T^,T]  and  all  vector  functions  f(t)  piecewise  continuous 

on  [T^,T],  the  following  Inequality  holds: 

T T 

71(T^,T)  - J f’(t)[J  K(t>T)f(T)dT]dt>0  0 


Lemma  [ 18]  ; For  the  class  of  kernels  K(t)  for  which  the  Laplace  transform 
exists,  the  necessary  and  sufficient  conditions,  in  order  that  K(t)  be  a 
positive  definite  kernel,  is  that  its  Laplace  transform  be  a positive  real 
transfer  matrix. 

Therefore,  by  the  above,  L,M,N,  and  S may  be  generalized  to  be 
positive  real  matrix  operators  Instead  of  constants. 
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APPENDIX  B 


The  following  assertion  provides  sufficient  conditions  under  which 

control  u will  achieve  perfect  model  following. 

P 

Assertion : Sufficient  conditions  for  the  control  (6)  to  force  the  plant  (4) 

to  follow  the  model  (5)  are: 

1.  (I  ■ • V ■ ° 

(I  - B b'*')B  “ 0 

ppm 

(I  - B b'*')R  ” 0 ; (I  - B b‘'’)R  - 0;  y-ltoP,  0=1  to  3 
ppm  P P P 


2.  lim  e - 0 
t “ 

A + A T -1  T 

where  e = X - X^  and  B ■ (B  B ) B 

m p p p p p 

Furthermore  a set  of  controls  that  will  achieve  perfect  following  under  the 
above  conditions  are  characterized  by: 


K + AK  - K = B (A  - A. ) 
m p p p m 1 


K + Ak  “ B B 
u u p m 


\ ■ S.  + % ■ “p 


m m 
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Proof:*  Note  that  perfect  model  following  implies  e*0  and  e“0,  and 

» 

conversely.  Therefore  we  will  assume  e*0  and  then  show  under  what 
conditions  e will  equal  zero.  To  this  end  we  define 


e - X -X 
ra  p 


and  subtract  (4)  from  (5)  with  the  control  (3)  with  P*1  and  P*1 


e=AX  +RF  +Bu  - AX  - RF  - B[(AK  - K)X 
mm  mm  mm  pp  pp  p p P P 


+ )F  + (K^  t )F^  + <K^  -F 


(B3) 


m m 


requiring  e=0  and  e=0  implies  x = x , therefore  (B3)  becomes 

m p 


e = 0 = [A  - A - B (K  + AK  - K )]X  + [R  - B (IC  + AlC  )]F 
m p p'  m P P P m p R R m 


+ I-  Bp  - Bp<Bj  + % )lFp  1B„  - B (K^  + BK„)lu 


tn 


(B4) 


P P 

but  sufficient  conditions  for  (B4)  to  hold  are  that  each  term  equal  zero  for 


all  X , F , F and  u , therefore 
p m p m 


A - A - B(K  +AK  - K)»0;  R - B(K„  +AK)*0; 
m p p'm  p p'^  m p'r  R'^ 


P P 
0 


m 


; B - B (K  + Ak  ) * 0 ; 

m p u u'^ 


- Bp  - B (Kp  + ) 

P P 

define  a set  of  perfect  following  controls.  Solving  for  K's 


(B5) 


K + AK  - K - b'*‘(A  -a  ) ; 1C  + Ak,  ■ b"^  R ; 

m ppp  nip  R R pm 


Kp  -f  4Kp  - -B-  R 


+ 

P ■ P 


K + Ak  ■ b b 
u u p m 


(B6) 


Uses  the  techniques  of  [29]  and  [25]. 


5 

! i 
i . 


! i 


Note  however  that  to  justify  the  use  of  the  pseudo-inverse,  we  have  to 
formulate  the  conditions  under  which  (B6)  satisfy  (B5)  [29] . Therefore 
substituting  (B6)  into  (B5)  we  get 
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« - “ i « ■ Vp>  ^ ^ ‘ 


(I  - B B ) R - 0 
P P'  P 


(I  - B B ) B - 0 ; 

ppm 


(B7) 


Note  also  that  when  P > 1 and/or  P > 1 we  get  similar  terms  to  the  above 
because  the  functions  distribute,  i.e. 

(I  - B B )■*■  R “0  ; (I  - B„B‘f)R  “ 0 ; 


P P 


P P 


V ■^''v  *<v  • \ -K\  • 

' Pof  ‘a 


m 


m 


y 


a ■ l,...,p  and  y - l,.,.,p 


Thus  we  have  shown  that  when  (B7)  is  satisfied,  (B6)  may  be  used  to  make 
e*0.  Therefore,  since  e is  assumed  to  equal  zero  and  e*0,  the  assertion  is 


established. 
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APPENDIX  C 
SIMULATION  PROGRAMS 

All  the  simulations  of  Chapter  6 were  made  on  a DEC-10  digital 
computer  using  the  C.S.M.P.  simulation  program.  This  program  is  quite  easy 
to  use  since  it  treats  the  elements  of  a system  block  diagram  as  a 
topological  network.  This  means  that  each  statement  is  given  a function 
type  (integration,  summer,  multiplier,  etc.)  and  the  statement  numbers  of  the 
operands.  Initial  conditions  and  parameters  for  the  functions  are  given 
after  all  function  statements  have  been  declared.  These  are  then  coded  by 
reusing  the  statement  number  of  the  particular  block  and  then  specifying 
parameter  values  to  the  inputs  of  the  blocks  by  a one-to-one  correspondence. 

The  two  programs  in  this  section  are  typical  of  all  the  simulation 
programs  made.  NLMR0.DAT,  for  example,  was  used  for  the  simulation  of  the 
error  signals  of  the  nonadaptive  system.  NLMRF.DAT  was  used  for  simulating 
the  adaptive  system  with  initial  conditions  and  adaptive  gains  of  10. 
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Table  2.  Program  NLMR0.DAT  Simulates  Nonadaptive  System 
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Table  3.  Program  NLMRF.DAT  Simulates  Adaptive  System  with  Initial  Conditions 
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